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We search for coincident gravitational wave signals from inspiralling neutron star binaries using LIGO and 
TAMA300 data taken during early 2003. Using a simple trigger exchange method, we perform an inter¬ 
collaboration coincidence search during times when TAMA300 and only one of the LIGO sites were opera¬ 
tional. We find no evidence of any gravitational wave signals. We place an observational upper limit on the rate 
of binary neutron star coalescence with component masses between 1 and 3 Mq of 49 per year per Milky Way 
equivalent galaxy at a 90% confidence level. The methods developed during this search will find application in 
future network inspiral analyses. 
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The first generation of gravitational wave interferometric 
detectors are rapidly approaching their design sensitivities. 
These include the LIGO Jj! and TAMA300 yj] detectors as 
well as GEO H and Virgo 01. Inspiralling binaries of neu¬ 
tron stars and/or black holes are one of the most promis¬ 
ing sources of gravitational radiation for these detectors. In¬ 
deed, several searches for such signals have already been com¬ 
pleted 0 6 , 0 , 8 ]. In the long term, the chances of detect¬ 
ing gravitational waves from a binary inspiral are greatly im¬ 
proved by making optimal use of data from all available de¬ 
tectors. The immediate benefit of a multi-detector coincidence 
search is a significant reduction in the the false alarm rate for 
a fixed detection efficiency. Additionally, a search involving 
all available detectors will provide an increase in observation 
time when, for example, at least two detectors are operating. 
The different orientations of the detectors make them sensitive 
to different parts of the sky, thus a combined search can lead 
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to improved sky coverage. If an event is detected in multiple 
instruments it is possible to localize the position of the source 
and improve parameter estimation. In addition, independent 
observations in well-separated detectors using different hard¬ 
ware and analysis algorithms would increase confidence in a 
detection, while reducing the possibility of an error or bias. 

The importance of joint searches has long been acknowl¬ 
edged, and indeed several network searches have previously 
been completed. A network of resonant decectors was used 
to carry out a joint search for gravitational wave bursts 0, 
and more recently data from the LIGO and TAMA300 detec¬ 
tors were used to perform a joint burst search O- In this 
paper, we present the first inter-collaboration search for grav¬ 
itational waves from the binary inspiral of neutron stars using 
modem large scale interferometric detectors. This represents 
an important step towards a global network analysis of grav¬ 
itational wave data. Furthermore, this search provides a firm 
basis for development of network analysis techniques. 

The joint coincidence search described here uses data from 
the second LIGO science run (S2) which occurred at the same 
time as the eighth TAMA300 data taking run (DT8) in 2003. 
The LIGO S2 data have already been searched for gravita¬ 
tional waves from binary neutron stars 0 . That search used 
only data in which both of the LIGO sites were operational. 
In this paper, we report on a coincidence search using LIGO 
and TAMA300 data when only one LIGO site was operating 
in coincidence with the TAMA300 detector. The LIGO data 
analyzed in this paper was not analyzed in Ref. 0. During 
S2 and DT8, the LIGO detectors were an order of magnitude 
more sensitive than TAMA300. However, since TAMA300 
was sensitive to the majority of candidate sources in the Milky 
Way, a joint coincidence search provides information about 
inspiraling neutron star binaries in the galaxy. Further, by 
performing this joint search between the LIGO and TAMA 
collaborations, we are able to significantly increase in the 
length of time searched in coincidence during the S2/DT8 run. 
Since LIGO and TAMA300 were the only large interferome¬ 
ters which were operated during S2/DT8 period, it is impor¬ 
tant to perform a joint analysis. 

The data from each of the detectors are searched indepen¬ 
dently for event candidates, or “triggers” 0. The details of 
these triggers, such as the coalescence time and the masses of 
the component stars, are then exchanged between collabora¬ 
tion members, and the triggers are searched for coincidences. 
The coincidence requirements of the search are determined 
by adding simulated signals to the data streams of the detec¬ 
tors, and determining the accuracy with which various param¬ 
eters are recovered |JT|]. The exchange of single instrument 
triggers and subsequent coincidence analysis is quite simple 
and does not involve the exchange of large amounts of in¬ 
terferometer data. It provides a natural first step in an inter¬ 
collaboration analysis. If an interesting candidate event were 
found, it would then be followed up by an optimal, fully co¬ 
herent analysis of the data around the time of the candidate. In 
this joint LIGO-TAMA300 search, we find no evidence of any 
inspiral signals in the data and so we place an observational 
upper limit on the rate of binary neutron star coalescence in 
the Milky Way galaxy. 


The LIGO network of detectors consists of a 4km interfer¬ 
ometer “LI” in Livingston, LA and a 4km “HI” and a 2km 
“H2” interferometer which share a common vacuum system 
in Hanford, WA. TAMA300 is a 300m interferometer “Tl” 
in Mitaka, Tokyo. Basic information on the position and ori¬ 
entation of these detectors and detailed descriptions of their 
operation can be found in Refs. il. The data analyzed in 
this search was taken during LIGO S2, TAMA300 DT8 be¬ 
tween 16:00 UTC 14 February 2003 and 16:00 UTC 14 April 
2003. We only analyze data from the periods when both LIGO 
and TAMA300 interferometers were operating. Furthermore, 
we restrict to times when only one of the LIGO sites was op¬ 
erational. Therefore, we have four independent data sets to 
analyze: the data set during which neither HI nor H2 were 
operating — the nHl-nH2-Ll-Tl coincident data set (here 
“n” stands for “not operating”) — and three data sets when 
one or both of the Hanford detectors were operational but LI 
was not — the Hl-H2-nLl-Tl, Hl-nH2-nLl-Tl, and nHl- 
H2-nLl-Tl coincident data sets. 

During the S2 science run, a strong correlation was found 
in the LI interferometer between inspiral triggers and non¬ 
stationary noise in the auxiliary channel, Ll:LSC-POB_I, 
which is proportional to the length fluctuations of the power 
recycling cavity. Therefore, we apply a veto to exclude times 
of excess noise in POB J, details of which are given in Ref. 
0. No efficient veto channels were found for the HI, H2 or 
Tl detectors. After applying the veto to LI, there are 34 hours 
of nHl-nH2-Ll-Tl data. Additionally, there are 334 hours of 
Hl-H2-nLl-Tl data, 212 hours of Hl-nH2-nLl-Tl data and 
68 hours of nHl-H2-nLl-Tl data, giving a total observation 
time of 648 hours. The data used in this search are summa¬ 
rized in Figure^ 

To avoid any bias from tuning our pipeline using the same 
data from which we derive our upper limits, the tuning of anal¬ 
ysis parameters was done without examining the full coinci¬ 
dent trigger sets. Instead, parameter tuning was done on the 
playground data which consists of approximately 10% of the 
data chosen as a representative sample. In this analysis, the 
length of playground data is 64 hours. The analysis of the 
playground data and tuning of the search is described in more 
detail in Ref. El- The playground data was searched for 
candidate gravitational wave events, but was excluded from 
the data set used to place the upper limit. Subtracting the play¬ 
ground data leaves a total of 584 hours of non-playground data 
used in placing the upper limit. 

In a search for inspiralling neutron star binaries, we can 
characterize the sensitivity of the detectors by their maximum 
observable effective distance, or range. This is defined as the 
distance at which an inspiral of 1.4-1.4 Mq neutron stars, in 
the optimal direction and orientation with respect to each de¬ 
tector, would produce a signal to noise ratio (SNR) of 8. The 
effective distance of a signal is always greater than or equal 
to the actual distance. On average it is 2.3 times as large as 
the actual distance, with the exact factor dependent upon the 
source location and orientation relative to the detector. During 
the S2 science run the ranges of the LIGO detectors, averaged 
over the course of the run, were 2.0, 0.9 and 0.6 Mpc for LI, 
HI and H2 respectively. This made them sensitive to signals 
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FIG. 1: The number of hours that each combination of detectors was 
searched during the S2/DT8 run. The upper number gives the amount 
of time the specific LIGO detectors were coincidentally operational. 
The lower number gives the total amount of time searched in coin¬ 
cidence with TAMA300. The shaded region corresponds to the data 
used in this search. 


from the Milky Way and favorably oriented potential sources 
in the local group of galaxies. The range of TAMA300 dur¬ 
ing DT8 was 52 kpc, making it sensitive to the majority of 
the Milky Way. Thus, the detectors were sensitive to a simi¬ 
lar population of candidate sources. Since we require a signal 
to be observed in both the LIGO and TAMA300 detectors, for 
this search we restrict our attention to gravitational waves pro¬ 
duced by inspiralling neutron star binaries in the Milky Way. 

The search methods employed in this paper are similar to 
those used in the LIGO S2 search 0 and the independent 
TAMA300 DT8 search 0 . Therefore, in this paper we 
will not describe the LIGO or TAMA300 analysis pipelines 
in great detail, but instead emphasize the differences between 
this search and those described previously. 

For the LIGO search, we split the data into analysis blocks 
of 2048 seconds length, overlapped by 128 seconds. For each 
block, we construct a template bank with a minimal match 
of 97% and component masses between 1 and 3 Mq 0. 
We analyze the data using the FINDCHIRP implementation of 
matched filtering for inspiral signals in the LIGO Algorithm 
Library 00. The most important thresholds used in the 
LIGO search are given in Table U Most notably, we use an 
SNR threshold p* = 7 for matched filtering. Additionally, 
we perform a waveform consistency (x 2 ) test 0. For this, 
we require the power observed in the signal to be evenly dis¬ 
tributed between p frequency bands. The threshold is 

X 2 <(P + $ P 2 )C ■ ( 1 ) 

We use a higher threshold on SNR (7 rather than 6) and also a 
tighter x' 2 threshold (5 rather than 12.5 in the Hanford detec¬ 
tors) than in the LIGO only S2 inspiral analysis. This is due to 
the fact that we limit our attention to signals from the Milky 
Way which tend to have a large SNR in the LIGO S2 data 
stream. The tighter thresholds vastly reduce the false alarm 
rate while giving a negligible loss of detection efficiency. 

For times during which both the HI and H2 detectors were 
operational, we perform a triggered analysis of H2, as de¬ 


Parameter Description 

value 

MM 

Templatebank Minimal Match 

97% 

P * 

Matched Filter Threshold 

7.0 

P 

Number of x 2 bins 

15 

5 

X' 2 threshold parameter 

0.023 

e 

\ 2 threshold parameter 

5.0 

StHH 

H1/H2 Timing Coincidence 

1.0 ms 

SrriHH 

H1/H2 Mass Coincidence 

0 


TABLEI: A list of the most significant parameters used for the search 
of the LIGO data. 


Parameter Description 

value 

MM 

Templatebank Minimal Match 

97% 

P* 

Matched Filter Threshold 

7.0 

P 

Number of X' 2 bins 

16 

8 

X 2 threshold parameter 

0.046 

r 

X 2 threshold paramater 

2.3 


TABLE II: A list of the most significant parameters used for the 
search of the TAMA300 data. 

scribed in detail in Ref. 0. We produce a template bank and 
matched filter the HI data. Only for those times and masses 
that we obtain a trigger in HI do we filter the H2 data. This 
significantly reduces our analysis time while having no effect 
on the detection efficiency. We then search for triggers co¬ 
incident in time and mass between the HI and H2 detectors. 
The use of a triggered search allows us to require the mass 
parameters of coincident triggers to be identical. Studies per¬ 
formed by injecting simulated signals show we can determine 
the end time of an inspiral to within 1ms and consequently 
we use this as our time coincidence window. Finally, we im¬ 
plement an amplitude consistency test between triggers in HI 
and H2 0; this includes keeping any triggers from HI whose 
recovered effective distance renders them unobservable in the 
less sensitive H2 detector. 

For the TAMA300 search, we split the data into analysis 
blocks of 52.4288 seconds length. The adjacent blocks of 
data are overlapped by 4.0 seconds in order not to lose signals 
which lie on the border of two adjacent blocks. We construct 
a template bank with a minimal-match of 97% III for each 
locked segment, in which the detector was continuously op¬ 
erated without any interruptions. The most significant thresh¬ 
olds in the TAMA300 search are listed in Table UTl We use a 
SNR threshold p* = 7 for matched filtering. In the TAMA300 
only search, we introduce a threshold on the value of p/ \fx 2 
to reduce the number of false alarms iTT^HTsll . However, in the 
LIGO-TAMA300 analysis, we introduce a x 2 threshold as in 
Eq. 0. By cutting on % 2 , the number of triggers is signifi¬ 
cantly reduced. In addition, some of the coincidence analysis 
becomes much simpler since LIGO and TAMA300 use a sim¬ 
ilar criterion for x_ 2 - More details of the TAMA300 analysis 
pipeline are available in Ref. 0111. 

The requirements for coincidence between triggers in the 
LIGO and TAMA300 detectors are determined by adding sim¬ 
ulated inspiral events to the data streams of the detectors. 
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Parameter Description value 

StHT Timing between Hanford and TAMA 27.0 ms 

5tLT Timing between Livingston and TAMA 35.0 ms 

5M Chirp mass window 0.05 Mq 

TABLE III: The coincidence windows used for the LIGO-TAMA300 
search. 


Thresholds are chosen so that injected signals seen separately 
in both the LIGO and TAMA300 detectors survive the co¬ 
incidence step with near 100 % efficiency, while minimizing 
the rate of accidental coincidences. Since both the LIGO and 
TAMA300 pipelines can accurately determine the coalescence 
time and mass of an injected signal, it is natural to require 
consistency of these values in our coincidence test. We mea¬ 
sure the accuracy with which these parameters are recovered 
in each detector and set the coincidence window to be the sum 
of these accuracies. The values of time and mass coincidence 
parameters are given in Table ED Both pipelines recover the 
end time with an accuracy of 1 ms, to which we must add the 
light travel time between sites to obtain the values given in the 
table. The mass parameter most accurately recovered by the 
pipelines is the chirp mass of a signal. The chirp mass is de¬ 
fined as A4 = Mr] 3 / 5 , where M = m\ + m 2 is the total mass 
of the system and 77 = is the dimensionless mass 

ratio. To pass coincidence, we require the chirp masses of two 
triggers to agree within 0.05 Mq. Further details of how these 
parameters were chosen are available in Ref. 0 . 

The coincidence parameters described above were chosen 
to provide a good efficiency to simulated events. However, 
there is some chance that noise induced events in the detec¬ 
tors might survive our coincidence tests. In order to estimate 
the background of such chance coincident triggers we per¬ 
form a time-shift analysis 0 - To do this, we time-shift the 
TAMA300 triggers by multiples of 5 seconds and search for 
coincidence between the time-shifted TAMA300 triggers and 
LIGO triggers. We perform 100 time-shifts, with a value of 
the time-shift ranging from —250 to 250 seconds. These shifts 
are much longer than the light travel time between the sites, 
so that any coincidence cannot be from actual gravitational 
waves. They are also longer than the typical detector noise 
auto-correlation time, longer than the longest signal template 
duration (4 seconds) and shorter than typical timescales of de¬ 
tectors’ non-stationarity, so that each time-shift provides an 
independent estimate of the accidental coincident rate. The 
SNRs of the triggers obtained from the time-shift analysis are 
plotted in Figure [2] The plot shows that the distribution of 
background coincidences does not follow the circular false 
alarm contours expected for Gaussian noise |20J. Instead, 
a statistic which more accurately reflects the constant false 
alarm probability contours is the sum of the SNR in the two 
detectors, 


Pc — Pligo T Ptama • (2) 

We use this statistic in our analysis to distinguish background 
triggers from detection candidates. 

To measure the sensitivity of the search, we perform a set 


co 

< 


14 

13 

12 



7 8 9 10 11 12 13 14 

LIGO SNR 


FIG. 2: The signal to noise ratios p U co vs Ptama of the accidental 
coincident triggers using 100 time-shifts. The contours of constant 
false alarm probability are also shown. 


of injections into both sets of data. The simulated waveforms 
added to the data consist of galactic binary neutron star inspi¬ 
ral signals. The majority of neutron stars in the Milky Way 
lie in the galactic bulge, which we take to have a radius of 4 
kpc and height of 1.5 kpc. The sun is assumed to lie 8.5 kpc 
from the center of the galaxy. Further details of the galac¬ 
tic model used are available in Ref. 12111 . The mass distri¬ 
bution is described in detail in Ref. |[22t l. Of the injections 
performed, 76% have an associated coincident trigger in the 
LIGO and TAMA300 detectors. The majority of the injections 
not detected have an effective distance at the TAMA300 site 
greater than TAMA300’s range during DT 8 . However, there 
were also a few injections which were very poorly oriented 
for the LIGO detectors, and hence have a large effective dis¬ 
tance, making them unobservable to LIGO. Finally, several in¬ 
jections produce triggers in both the LIGO and TAMA300 de¬ 
tectors but these fail our coincidence requirements. The SNRs 
of these triggers are close to threshold in TAMA300 and the 
injection parameters, in particular the chirp mass, are recov¬ 
ered poorly. In Figure^we plot the coincident triggers associ¬ 
ated with injections superimposed on those from the time-shift 
analysis. This shows that triggers from the found injections 
are well separated from the accidental coincidences found in 
the time-shift analysis. 

In Figure^] we plot the sensitivity of the search to injected 
Milky Way signals. For consistency with previous searches 
0 ] we use Nq to represent the number of galaxies to which 
the search is sensitive. For this search, Nq is equivalent to 
the fraction of Milky Way signals we are sensitive to. The 
figure shows the number of galaxies the search is sensitive to 
as a function of the threshold on the combined statistic given 
in Eq. Thus, at threshold we are sensitive to a little more 
than three quarters of candidate sources in the galaxy. The 
efficiency curve is used later in determining the upper limit 
and associated systematic errors. 
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FIG. 3: The signal to noise ratios p L ioo vs Ptama of the triggers as¬ 
sociated with injections (+) and those from accidental coincidences 
arising in 100 time-shifts (o). 



FIG. 4: The efficiency of the LIGO-TAMA300 joint analysis to 
simulated galactic inspiral events. The number of galaxies ( Ng ) to 
which the search is sensitive is plotted as a function of the threshold 
on the combined statistic p c (= puoo + Ptama). 


We analyze the S2/DT8 data using the pipeline described. 
The cumulative distribution of p c of the coincident triggers 
is shown in Figure 0 On this plot, the expected number of 
triggers obtained from the time-shift analysis is shown, as well 
as the standard deviation of the number of triggers obtained in 
the time-shifts. The results of the analysis of the full data 
are overlayed on top of this. It is clear from the figure that 
the distribution of coincident triggers is consistent with the 
background estimated from time-shifts. There are no triggers 
with combined SNR greater than p max = 15.3. Therefore, 
we conclude that there is no evidence for gravitational wave 
signals in the LIGO-TAMA300 S2/DT8 data set. 

Given the set of triggers displayed in Figure|3we can obtain 


FIG. 5: The triggers from the analysis of the full LIGO-TAMA300 
data set. The x represent the expected background number of trig¬ 
gers at or above a given combined SNR p c based on the 100 time- 
shifts performed. The bars indicate the standard deviation of the 
number of events, calculated from the time-shift results. The trig¬ 
gers from the final S2/DT8 data set are shown as o. 


an upper limit on the rate of binary neutron star coalescences 
per year per Milky Way Equivalent Galaxy (MWEG). (Al¬ 
though this search is only sensitive to galactic inspiral events, 
we maintain the standard “MWEG” H for describing the 
upper limit). We use the loudest event statistic |2,,3|. which 
makes use of the detection efficiency at the combined SNR of 
the loudest event in order to construct the upper limit. The 
90% confidence frequentist upper limit is given by 


_ 2.303 +In P b 
90% “ T7V G (p max ) ' 


(3) 


In the above, T is the observation time of 584 hours, Pb is the 
probability that all background triggers have a SNR less than 
Pmax? and Nq is the number of MWEGs the search is sensitive 
to at the combined SNR of the loudest event p max . Nq is de¬ 
termined from Figure 0] to be 0.76 MWEG for p max = 15.3. 
Although the time-shift analysis provides us with an estimate 
of Pb = 0.2, we note that it is difficult to establish a system¬ 
atic error associated with this estimate, and therefore take the 
conservative choice of setting Pf, = 1. From these numbers, 
we obtain an upper limit of P go % = 45 y -1 MWEG^ 1 . 

The possible systematics which arise in a search for bi¬ 
nary neutron stars are described in some detail in Ref. 0, 
and we will follow the analysis presented there to calculate 
the systematic errors for the above result. The most signif¬ 
icant effects are due to the possible calibration inaccuracies 
of the detectors, the finite number of Monte Carlo injections 
performed, and the mismatch between our search templates 
and the actual waveform. We must also evaluate the system¬ 
atic errors associated with the chosen astrophysical model of 
potential sources within the galaxy. All systematic effects in 
the analysis pipeline (such as less than perfect coverage of the 
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template bank) are taken into account in the Monte Carlo es¬ 
timation of the detection efficiency. 

This search was sensitive to most, but not all, signals from 
the Milky Way. Thus, the specific model of the source dis¬ 
tribution within the galaxy will affect the upper limit. The 
majority of the mass in the galaxy, and hence the potential 
sources, is concentrated near the galactic center. Therefore, 
our efficiency will be most affected by changing the distance 
from the sun to the center of the galaxy in the model. In 
this search, the sun’s galactocentric distance is assumed to be 
8.5 kpc. Varying this distance between 7 and lOkpc leads 
to a change in efficiency of 0.04 MWEG. Different mod¬ 
els for NS-NS formation can lead to variations in the NS 
mass distribution. Based on simulations with a 50% reduc¬ 
tion in the number of binary systems with masses in the range 
1.5M 0 < 7711,7712 < 3.0Mq, we can estimate the variation 
in Nq to be 0.01 MWEG 

Any calibration inaccuracy in TAMA300 could have a sig¬ 
nificant effect upon our efficiency. This is clear from Fig- 
ure 0 which shows a significant number of injections found in 
TAMA300 close to threshold. Two effects contribute to this 
calibration error: an overall normalization error (associated 
with the magnetic actuation strength uncertainty and its effect 
on calibration), and uncertainty in the frequency-dependent 
response. The error in the normalization is of order 5%, but 
the long-term drift is unknown, so we conservatively use 10% 
in this paper. The frequency-dependent error was estimated 
and shown to be <C 10%, so it is subsumed into the overall 
10% error on the SNR of the triggers. This calibration un¬ 
certainty leads to a 0.02 MWEG effect on our efficiency. The 
majority of injections are observed well above threshold in the 
LIGO detectors, and consequently the calibration uncertainty 
of 8.5% in LI and 4.5% in H1/H2 results in a smaller uncer¬ 
tainty in the efficiency of < 0.01 MWEG. The error in the 
efficiency measurement due to the finite number of injections 
performed is 0.01 MWEG. Differences between the theoreti¬ 
cal waveforms used in matched filtering the data and the real 
waveforms would decrease the efficiency of our search. Al¬ 
lowing for a 10% loss in SNR due to inaccuracies in the model 
waveform HIS El leads to a +0/ - 0.02 MWEG effect 
on the efficiency. Combining these effects, we obtain an effi¬ 
ciency of Nq = 0.76^o Qg. Taking the downward excursion 
on Nq, we obtain a conservative upper limit of 

7 Z 90% = 49 y -1 MWEG -1 . (4) 

This rate is substantially higher than the predicted astro- 
physical rate of 8.3 x 10 -5 y -1 MWEG -1 |27j. However, 
the rate limit obtained in this paper is comparable with the 
rate limit of 47 y -1 MWEG -1 obtained from the LIGO-only 
S2 search 0 , which was performed on a complementary data 
set. Since these searches were performed on independent data 
sets, if astrophysically relevant, these upper limits could then 


be combined to produce the best possible limit. The fact 
that the LIGO S2 and LIGO-TAMA300 S2/DT8 limits are so 
similar demonstrates that the overall sensitivities of the two 
searches are very nearly equal. This is achieved despite the 
fact that the TAMA300 detector was less sensitive than LIGO 
during S2/DT8. The high duty cycle of TAMA300 (over 80%) 
compensates for the reduced sensitivity, and leads to a similar 
overall result. 

In this paper, we have presented the methods and results 
from the first multi-collaboration, network search for gravi¬ 
tational waves from inspiralling binary systems using large 
scale interferometers. The search was performed using a trig¬ 
ger exchange method, requiring coincidence in both the end 
time and chirp mass of triggers between instruments. Us¬ 
ing this method, we have performed all necessary steps of 
the analysis, including time-shifts, signal injections and the 
calculation of the upper limit. The joint, coincidence search 
presented here is a natural first step in any network analysis. 
The methods developed during this search will be applied in 
future network searches. Indeed, the experience gained during 
this joint search is being used in subsequent LIGO Scientific 
Collaboration searches of LIGO and GEO data. Furthermore, 
a trigger exchange coincidence analysis is being developed 
as the first stage of a future joint LIGO-Virgo analysis. The 
optimal network search would likely involve a fully coherent 
analysis El of the detectors’ data streams around the times 
of coincident events. A coherent followup to the coincidence 
method presented in this paper would be included in future 
analyses of LIGO and TAMA300 data with improved sensi¬ 
tivity, or in joint analyses of the planned second generation 
detectors such as advanced LIGO 128] in U.S. and LCGT ll29tl 
in lapan. 
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